Non-degenerate four-wave mixing in rubidium vapor: transient regime 
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We investigate the transient response of the generated light from Four-Wave Mixing (FWM) in 
the diamond configuration using a step-down field excitation. The transients show fast decay times 
and oscillations that depend on the detunings and intensities of the fields. A simplified model taking 
into account the thermal motion of the atoms, propagation, absorption and dispersion effects shows 
qualitative agreement with the experimental observations with the energy levels in rubidium (5Si/2, 
5-Pi/2; 5P3/2 and 6S1/2). The atomic polarization comes from all the contributions of different 
velocity classes of atoms in the ensemble modifying dramatically the total transient behavior of the 
light from FWM. 

PACS numbers: 42.65.-k 32.80.Qk 



I. INTRODUCTION 

Four- Wave Mixing (FWM) is a non-linear optical pro- 
cess that coherently combines three waves to produce a 
fourth one, and is used in many areas of physics. In 
atomic physics FWM is utilized to study spectroscopic 
and quantum optics processes [lj-Q . It can generate spe- 
cific frequencies not available with direct electromagnetic 
sources, and has become a widely used tool for condensed 
matter physics to characterize the short time correla- 
tions and coherences of semiconductors (see for exam- 
ple Refs. 0, 0, 11). Recent work with FWM for quan- 
tum information includes implementations of the atomic- 
ensemble-based quantum repeater of Duan, Lukin, Cirac, 
and Zoller M such as |10l - fl2j , the generation of strong 
squeezing [l3|, [l4[ and quantum imaging [15j . These 
implementations rely on the coherences created in the 
atomic systems, which are affected by dephasing mech- 
anisms producing decoherence over time. Study of the 
temporal dynamics of FWM can provide insight into the 
decay of the atomic polarization and the process of de- 
coherence. 

We study the time evolution of the atomic coherences 
in FWM in rubidium vapor with transitions at telecom- 
munication wavelengths. We use a step-down excitation 
of the excited states and observe the generated light from 
FWM in a diamond configuration, where the levels in- 
clude the 5S ground state, the 5P states (Dl and D2 
lines), and the 6S state ( See Fig. [T]). This study com- 
plements the investigation of the frequency response of 
the process in Refs. [H, [13] and concerns effects oth- 
erwise unobservable in the frequency regime due to the 
nonlinearity of the system. The transient signal contains 
information about the temporal evolution of the induced 
atomic polarization in the medium and can have collec- 
tive atomic effects [HI, [l9| . Since the atomic polarization 
is responsible for quantum effects understanding the 
decoherence is fundamental for app lications of FWM in 
quantum information; our work [2CJ shows the limitations 
on non-classicality from decoherence. 



The FWM signal in our configuration in steady state 
is weak. We step down to zero the laser intensity that 
connects one of the intermediate states to the highest ex- 
cited state and observe the time evolution of the FWM 
as it decays and propagates out of the system. The time 
scales of this atomic system allow us good resolution with 
this simple step-down process; this is in contrast to the 
situation encountered in FWM in semiconductors where 
the dynamics are very fast (picoseconds or faster) and 
pump-probe techniques such as time-integrated FWM 
and time-resolved FWM are used to study correlation 
dynamics in semiconductors [2l| - |25| . 

The organization of this work is as follows: Sec. |TT] 
describes a simplified theoretical model of the system. 
Sec. Mil contains the experimental setup and observa- 
tions. Sec. IIVI contains the analysis incorporating other 
effects present in the experiment such as Doppler broad- 
ening, and Sec. |V]explains the results of the investigation 
of the transients in FWM. Sec. I VII summarizes our con- 
clusions. 



II. THEORY 

We study the evolution of the atomic coherences cre- 
ated by the non-linear interaction of three fields in the 
atomic medium in the FWM process when one of those 
fields is suddenly set to zero. The density matrix treat- 
ment allows us to study the decay and possible dephasing 
mechanisms such as Doppler broadening. 

Figure Q] shows the fields Ei at frequencies u>i coupling 
the four levels of the atomic system. E\ couples the tran- 
sition I a) — > I b), E2 the transition \b) — s> |c) and £3 the 
transition \a) — > \d). The generated light from FWM is 
resonant with the transition \d) — > |c). We assume that 
the generated light from the FWM process is weak in 
comparison with the incident fields and we ignore prop- 
agation effects at this frequency. 

The density matrix equations describing this system in 
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FIG. 1: (Color online) Four-level atom interacting with three 
external fields Ei, E2 and E3 with frequencies wi (795 nm), C02 
(1324 nm) and ^3 (780 nm) respectively, and generated light 
(red online) via four- wave mixing at frequency 0J4 (1367 nm). 
Ei and E3 are continuous wave (CW) fields and E2 is time 
dependent. The energy levels in the diamond configuration 
in 85 Rb are \a) = 5S 1/2 , F = 3, [6) = 5P 1/2 , F = 2, |c) = 
6S 1/2 , F = 3 and |d) = 5P 3/2 , F = 4. 



the rotating wave approximation (RWA) [4| are: 

Pfi = ~ ^ r ni/ 9 n + ^ ^nlPnn (1) 

n n 

— r ^ (VJmPmi — P/m^mi) 
m 

Pmn — {j'^mn Tmn) Pmn 

V m ) 

where pu describes population of level I, p mn the coher- 
ence between the levels m and n, T mn the decay rate 
from level m to n, T n = l/r n the total decay rate of level 
n with r n the lifetime of this level and j mn = |(r m +r„) 
the dephasing rate of the coherence p m n- The term A mn 
is the detuning of the external field coupling the transi- 
tion \n) — > \m) from the energy of this transition (see Fig. 
[1]). Energy conservation in FWM, cji +CJ2 = ^3+1^4, de- 
fines the frequency and detuning for the generated light. 
The term V nm = —Ml nm describes the interaction of the 
atomic ensemble with the external field in the dipole ap- 
proximation U where O nm is the Rabi frequency for the 
external field coupling to the transition \n) — > \m) given 
by fli — pi ■ Ei/h, where pi is the electric dipole moment 
of the transition \n) — > \m) and Ei is the electric field 
coupling to this transition. 

The density matrix equations are 16 coupled equations. 
Conservation of population, Pnn — 1, reduces this sys- 
tem to a set of 15 coupled differential equations. The 
generated light from FWM resonant with the transition 
\d) — > |c) (see Fig. [T]) is proportional to the atomic po- 
larization oscillating at the frequency W4 (l6| . The in- 
duced atomic polarization for the transition \d) —y |c) 



is P<fc = NpdcPdc where N is the atomic number den- 
sity 0, HH . The study of the transient behavior of the 
light from FWM requires solving these coupled differen- 
tial equations for the atomic coherence pdc- 

We assume that the system starts in the steady state 
at time t = when we turn off the field E2 (at t = 
0). The system evolves in the presence of only E\ and 
E3 with the steady state as the initial conditions. This 
allows some simplifications in the calculation. We obtain 
a decoupled subset of the system of 15 coupled equations 
that contains the term pd c - 

r l I Pea \ I &ca 

ifll ifl 3 \ I p ca \ 

— Pcb = i^i a cb p cb (2) 

\Pdc I \ 1^3 Ot dc ) \ Pdc ) 

where a dc = -iA FWM - led, a ca = -i(Ai + A 2 ) - 7 CQ 
and a co = — iA 2 — "f c b with Aj the detuning of the laser 
from the corresponding atomic transition, as shown in 
Fig. [l] and Afwm = Ai + A 2 — A 3 . The steady state 
solution of the density matrix equations determines the 
initial phases and amplitudes for these coherences. The 
solution for pd c is a linear combination of three exponen- 
tial functions of time with different characteristic energies 
(quasienergies), 

Pdc = A ie -^ + A 2 e-^ t + A 3 e~ lX3t (3) 

where Aj (j — 1,2,3) are complex time-independent am- 
plitudes and the quasienergies Aj are complex quantities 
whose real parts describe the oscillatory behavior of the 
generated light and imaginary parts the exponential de- 
cay constants. These quantities depend on the detunings 
and the Rabi frequencies of the lasers, and the ampli- 
tudes Aj are determined from the initial (steady state) 
conditions of the system. 

The real parts of the quasienergies show avoided cross- 
ings due to the interaction of the external fields with the 
atom. Ref. [16( describes our investigations of the FWM 
resonant structure where we compare those quasiener- 
gies with the measured spectra. That work, using the 
Schrodinger equation approach, does not take into ac- 
count dissipation. With the density matrix formalism it 
is now possible to see the damping of the atomic polar- 
ization. The intensity of the generated light from FWM 
is proportional to |pdc| 2 - The quantity |pdc| 2 j using Eq. 
©, contains three squared amplitudes associated with 
the same quasienergy, 

F j = \A j \ 2 e- 2I ™W, (j = 1,2,3) (4) 

showing exponential decay with a decay constant of 
2Im[Aj], and three cross terms plus their complex con- 
jugates (c. c), associated with different quasienergies, 

Gj = A 3 A* k e- l[x >- x * ]t + c. c, (5) 

for (j, k = 0,1,2), j 7^ k. The latter show a damped 
oscillatory behavior with an oscillation frequency (3jk = 
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FIG. 2: (Color online) FWM transients from the decoupled 
model, a) Square amplitudes Fj from Eq. Q. b) Damped 
oscillations Gj from Eq. Dotted (red) lines correspond 

to the quasienergy Ai, dashed (black) lines to A2 and solid 
(blue) lines to A3, c) Transient signal from the decoupled 
model using Eq. @ (circles) and numerical simulation using 
the density matrix Eqs. {TJ (solid line). Parameters used in 
the calculations: (Ai,A 2 ,A 3 ) = (0,270,30) MHz and Rabi 
frequencies of (Oi, f2 2 , f2 3 ) = (50, 10, 10) MHz. 



Re[Aj — AJ] and an exponential damping constant Qk = 
Im[Aj — A£]. Figs. [2^l and[2p show the quantities Fj and 
Gj as a function of time. The calculations use the rele- 
vant parameters for the experiment. The relevant levels 



are \a) = \5S 1/2 ), \b) = |5P 1/2 ), \c) = 
|5P 3 / 2 }, with lifetimes r 5 p 1/2 = 27.8 ns, 



in 85 Rb (Fig. EJ 
I61S1/2) and \d) = 

T5P 3/2 = 26.4 ns and tqs 1/2 — 46 ns and with correspond- 
ing linewidths r 5Pl/2 = 5.8 MHz, r 5 p 3/2 = 6 MHz and 
r 6 s 1/2 = 3.2 MHz. Figured: shows the transient FWM 
light using the simplified model from Eq. ([3]) and that 
numerically calculated from the complete set of density 
matrix equations Eqs. (HJ. 

For this particular set of Rabi frequencies and detun- 
ings, the transient signal of the generated field shows 
an oscillatory exponential decay behavior. The oscilla- 
tion frequency depends on the strength of the oscillatory 
terms from Eq. ([5]) . The decay of the total signal is domi- 
nated by the strongest components of the decaying terms 
in Eq. ((4]) . The shape of the transient signal depends on 
these six terms with the phases of the oscillations deter- 
mined by the initial conditions. The oscillation frequen- 
cies depend on the detunings and the Rabi frequencies of 



the fields. 

The decoupled subset of equations describes the be- 
havior very well, as seen in Fig. with the difference 
between the two calculations of the order of one part in 
10 3 , which is within the numerical error. This allows 
us to study the physical mechanism acting in this sys- 
tem using the simpler set of equations. The atomic co- 
herences in the decoupled system can be represented as 
three coupled oscillators. The oscillators corresponding 
to the atomic coherences pdc between c and d and p c b be- 
tween c and b can store and exchange energy because of 
their mutual coupling. This exchange process depends on 
their relative phases. The coherence p ca between c and 
a acts as part of the coupling coefficient in the system of 
two coupled oscillators pdc and p c b- Thus it produces a 
dynamic coupling between pd c and p c b- 



III. EXPERIMENTAL OBSERVATIONS OF 
FWM TRANSIENTS 



The specific hyperfine levels that the fields couple in 
the experiment are such that a field at 795 nm with fre- 
quency CJ795 (wi) couples the transition 55i/2, F = 3 — >• 
5-Pi/2> F — 2, a field at 1324 nm with frequency 
W1324 (102) the transition 5Pi/2, F = 2 — > 6Sx/2, F = 3 
and a field at 780 nm with frequency W780 (^3) the tran- 
sition 5S 1/2 , F = 3 ->■ 5P 3/2 , F = 4 in 85 Rb (see Fig. EE). 
The generated light at 1367 nm from FWM is resonant 
with the transition 5P 3 / 2 , F = 4 -4 65*1/2, F = 3 in 85 Rb 
with a frequency CJ1367 (0J4). The estimated Rabi frequen- 
cies of the external excitations (40% uncertainty due to 
their intensity profile) are ^795 = 45 MHz, 5I1324 = 8 
MHz and fl 78 o = 8 MHz. 

We use a fast fiber-coupled electro-optical modulator 
(EOM) to extinguish the excitation at 1324 nm in less 
than 4 ns, which is shorter than any atomic lifetime of 
the levels involved in the FWM process. A single photon 
detector collects the light from FWM as we turn off the 
1324 nm light while the 780 nm and 795 nm light beams 
are operating in continuous wave (CW) mode. 

Figure [3] shows the experimental configuration for the 
study of the transient regime of the FWM after the ex- 
tinction of excitation light at 1324 nm. This setup is 
similar to that used in the study of the process in the 
frequency domain described in Ref. [l6| , but here we use 
a rubidium cell at 100°C containing enriched 85 Rb (99% 
purity specified by the manufacturer) without buffer gas 
and without coating on the windows. The cell is L =1.5 
cm long and 2.5 cm in diameter. The typical powers of 
the lasers are: P 795 =2 mW, P 780 =20 ^W and the 1324 
nm has a power of Pi3 2 4=0.8 mW in the CW mode af- 
ter the fiber modulator. All the beams have horizontal 
polarization as does the generated light from FWM. The 
angles chosen to satisfy the phase-matc hing condition for 
the FWM process are shown in Fig. [3] [16j. We stabilize 
the frequencies of the lasers at 780 nm and 795 nm using 
RF sidebands and saturated absorption to lock their fre- 
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FIG. 3: (Color online) Schematic of the experimental con- 
figuration for the study of transient FWM. An electro-optic 
modulator (EOM) produces a pulse of light at 1324 nm with 
a turn-off of less than 4 ns while two lasers at 780 nm and 
795 nm work in continuous wave (CW) mode. A triggered 
single photon detector detects FWM light as a function of 
delay from the excitation pulse. The angles among the beams 
satisfy phase matching for the FWM. 



quencies using the absorption lines in rubidium. We tune 
the frequencies of the lasers using acousto-optic modula- 
tors in the frequency locking setups. We use a transfer 
cavity lock [27j for the frequency stabilization of the laser 
at 1324 nm. 

A 10 Gbits/s Lithium Niobate Electro-Optic Modula- 
tor (EOM) from Alcatel-Lucent with a bandwidth of 10 
GHz modulates the amplitude of the 1324 nm light. A 
digital delay generator SRS DG535 produces a 300 ns 
pulse of 4 Volts and provides the AC signal driving the 
EOM. We use a DC power supply for the DC voltage to 
the EOM for maximal suppression of light in the off stage 
of about 98% ± 2% at 1324 nm. 

The Single Photon Benchtop Receiver from Princeton 
Lightwave, model PGA-600, is an In-GaAs based single 
photon avalanche photodiode (APD) running in a gated 
mode to detect the generated FWM light. The gate has 
a 1 ns width which defines the temporal resolution in the 
experiment. We use optical attenuators to attenuate the 
light from FWM going into the fiber-coupled detector to 
keep a maximal rate less than 1 x 10 -2 photons per 1 ns 
window. We use a FEL 1350 long-pass filter from Thor- 
labs before a single mode fiber (SMF) with transmission 
of 70% at 1367 nm and no appreciable transmission at 
1324 nm. (We do not observe any counts caused by the 
1324 nm laser when we trigger this detector at a rate of 
100 x 10 3 counts per second). We use Lab View program- 
ming and GPIB and serial communication to control the 
SRS DG535 and the PGA-600 for the generation of 1324 
nm optical pulses and data collection of the light from 
FWM. 

We study the transients in the light from FWM as a 
function of powers of the CW lasers and detunings of the 
excitation fields from specific atomic transitions in 85 Rb 
in a diamond configuration as shown in Fig. [TJ The 
detunings and powers determine the exact values of the 
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FIG. 4: (Color online) a) Excitation pulse at 1324 nm (in 
black) and light from FWM at 1367 nm (in blue), b) Decay 
time of the transient FWM as a function of 780 nm power 
for temperatures T=344 K in (black) points and T=388 K in 
(red) triangles for P795 = 2 mW. b has its vertical axis zero 
suppressed. 



complex quasienergies that govern the atomic polariza- 
tion (Eq. ((3])). We observe pure decay for certain regions 
of power and detunings (see Fig 2J but also oscillatory 
decay for other regions (see Fig. [SJ. 

Figure 0Ji shows the excitation light at 1324 nm (in 
black) and the light from FWM at 1367 nm (in blue) as 
a function of time when the light at 1324 nm turns off at 
t = for the detunings (A 795 , Ai 324 , A 780 ) = (0, 100, 80) 
MHz and laser powers (P 795 , P 1324 , P 780 ) = (2,0.8,0.02) 
mW. Note the fast decay time of the light from FWM af- 
ter the extinction of light at 1324 nm. Fig. Hh> shows the 
decay time extracted from least squares fits to simple ex- 
ponential functions of the light from FWM as a function 
of the power of 780 nm light. The errors are calculated 
from the change of the extracted parameter to increase 
the minimized \ 2 of the fit by one. 

We observe that the decay time decreases as we in- 
crease the power of the 780 nm light, which might point 
to dephasing from power broadening of the level 5P3/2 
affecting the coherence pdc- This decay time is indepen- 
dent of atomic density in a range of 7.89 ± 2.51 x 10 11 to 
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1.5±.35x 10 13 atoms/cm 3 corresponding to temperatures 
T = 344 ± 5°K to 388 ± 5°K respectively as seen in Fig. 
HJd. This suggests that the short decay time is not due to 
any superradiant effects [28| or superfluorescence in the 
FWM process [19| . The short decay time of the transient 
FWM can be caused by dephasing mechanisms related to 
atomic diffusion and Doppler broadening in our sample 
[29j | . The decay does not depend on the power of the 795 
nm light in a range from 20 fiW to 2 mW. 

When the external fields interact with the atomic sys- 
tem, they generate an atomic polarization grating in the 
atomic ensemble which can generate light from FWM 
30]. This grating can diffuse due to the thermal motion 
of the atoms in the ensemble and contribute to dephas- 
ing in the process, shortening the decay time [20]. The 
study of this dephasing process as a function of thermal 
motion would require larger changes in the temperature 
of the ensemble than the temperature range used in this 
experiment. 

We study the transient behavior as a function of de- 
tunings of the three different lasers. We use 20 /iW of 
power for the 780 nm light and 2 mW of power for the 795 
nm light. The power of the 1324 laser in the CW mode, 
i.e., full transmission through the EOM, is 0.8 mW. This 
set of powers results in a FWM transient signal with a 
relatively long decay time, TExp — 6.0 ns, as shown in 
Fig. |4Jd, and the maximum generation of FWM light. 
This allows for a more detailed study of the transients in 
FWM as a function of detunings. 

Fig. [5] shows examples of the experimental observa- 
tion of the light from FWM in the transient regime as a 
function of the detuning of the 795 light, for A780 = 30 
MHz and Ai324=270 MHz. We observe different features 
in the generated light: a FWM-spike as we turn off the 
1324 nm for A 795 = -100 MHz and A 795 = -60 MHz, 
and oscillations at different frequencies. The FWM spike 
has an intensity above the steady state level, which corre- 
sponds to negative times in the figures. The decay times 
of the transients shown in Fig. [5] are approximately 6 ns 
which is consistent with the decay times in Fig. |4Jd for 
P780 = 20 fiW. The oscillation frequency in FWM de- 
pends on the atomic polarization ( See Eq. © and Fig. 
[5]) This is different from quantum beats from the hyper- 
fine structure [20j , whose frequency remains mostly fixed 
to power and detuning changes. 

This effect has been observed in the study of the ex- 
change excitation of a system of N two-level atoms cou- 
pled to a cavity mode in the linear [3l| and nonlinear (32| 
regimes, where oscillations and spikes in the transient 
regime are the result of exchanging energy between the 
electromagnetic field and induced atomic polarizations 
(vacuum Rabi oscillations [33|, HH). The asymmetry of 
the FWM-spike in our experiment comes from the phases 
in the atomic polarizations that depend on the detunings 
and the intensities (Rabi frequencies) of the three fields 
at 795, 780 an 1324 nm. 

We have explored the light from FWM for different 
regimes of laser detunings, and we have observed similar 




FIG. 5: (Color online) Experimental FWM transients as a 
function of 795 detuning for A780 = 30 MHz and Ai324=270 
MHz. a) A795 = -100 MHz. b) A795 = -60 MHz. c) 
A795 = MHz. d) A795 = 100 MHz. The excitation pulse 
is extinguished at t=0 after the system has reached steady 
state. Continuous lines (blue) join the data points (red). 



features in the transients as those shown as a function of 
795 detuning in Fig. [5] 

The parameter space in detuning and Rabi frequen- 
cies for FWM is large. Our exploration tries to highlight 
three main features: damping faster than the lifetime 
of the transitions (Fig. 0} when the detuning are large, 
while for smaller detuning we look into the overshoot 
from energy stored in the atomic polarization (Fig. [S]) 
and the oscillations from the atomic polarization inter- 
ference (Fig. [5J or a combination of both. At the end of 
section [V] we map out most of the region explored with 
respect to the predictions of the model of Eq. © . 



IV. ANALYSIS 

When we consider a thermal atomic gas, the atoms 
will have velocities described by the Maxwell-Boltzmann 
distribution (35[. Different velocity classes of atoms will 
observe different frequencies of the external fields due to 
the Doppler effect [35| . The transient response in FWM 
depends on the detunings of the lasers and also on the 
velocity of the atoms in the ensemble. 

The Doppler effect changes the transient behavior of 
the light dramatically when we consider the tempera- 
ture of the atomic ensemble. We take into account the 
Doppler broadening by introducing velocity-dependent 
detunings, and averaging the density matrix elements 
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FIG. 6: (Color online) Calculation of the FWM light in the 
transient regime after E2 is turned off at time t=0 as a func- 
tion of Ai (A795), for A2 = 270 MHz, corresponding to A1324, 
and A3 = 30 MHz, corresponding to A780, with Rabi frequen- 
cies (ni,n 2 ,n 3 ) = (95,20,20) MHz. The left column corre- 
sponds to the Doppler-free case (n=0); the right column after 
Doppler averaging. 



over the Maxwell-Boltzmann distribution [35[ 

Prnn = \ Pmn{v)F '(v)dv (6) 



where P(v) = v / m/2Trk B T exp [~mv 2 / (2k B T)} is the 
one-dimensional Maxwell-Boltzmann distribution with m 
the atomic mass, k B the Boltzmann constant and T the 
temperature of the gas. The total generated light from 
FWM is the result of contributions from all the different 
velocity classes in the ensemble. 

Figure [6] shows calculations of the transient behavior 
of the light from FWM for two detunings Ai (A795), 100 
MHz and -100 MHz, with A 2 = 270 MHz, and cor- 
responding to A1324, A 3 = 30 MHz, corresponding to 
A 780 , with Rabi frequencies (fii,0 2 ,^3) = (95,20,20) 
MHz. The left column corresponds to zero atomic ve- 
locity; the right column corresponds to the case includ- 
ing the Doppler effect due to the motion of the atoms 
at T =100°C using Eq. © and the 1-dimensional (1- 
D) Doppler average from Eq. (j6j). The different veloc- 
ity classes in the ensemble contribute to the light from 
FWM and modify its lineshape from the Doppler-free 
case. Note the shortening of the decay time of the tran- 
sients after the Doppler averaging. 

The oscillation frequencies of the transient FWM 
should depend on the detunings of the lasers as well as 
on the strength of the Rabi frequencies f2i and f2 3 . They 
do not depend on the Rabi frequency VL2 since it does not 
enter into the decoupled equations, Eqs. ([!]), because it 



is off; only contributes to the initial phases of the os- 
cillations. We use the steady state resonant structure of 
the FWM light from the experiment to estimate the Rabi 
frequency strengths. 

We next consider absorption of the light at 780 and 
795 nm from hyperfine levels in Rb, neglecting any ab- 
sorption at 1324 nm. The absorption modifies the Rabi 
frequencies entering into the process for different wave- 
lengths of the external fields. 

The linear absorption and dispersion of light in an 
atomic medium is described by the first order suscep- 
tibility x where = x^E. For a two-level atom 
1: 

X W(A)=A-1-, (7) 
A — 11 

where A is the detuning of the laser from the atomic 
transition, T is the linewidth of the transition and A is 
a constant. We assume that the total absorption of light 
from different hyperfine levels can be described by a to- 
tal linear susceptibility equal to the sum of the linear 
susceptibilities of each transition, 

where Ai defines the relative transition strengths of the 
transitions and A; the detuning of the laser from each 
hyperfine state. We account for the Doppler effect by 
averaging over velocities to obtain the total Dopplcr- 
averaged linear susceptibility XrltaV wnere the bar de- 
notes Doppler-averaged. The field propagating in the 
medium is d, H|| : 

E(z) = E(0) exp (-i0D(z) MXrltJ) (9) 

where E(0) is the incident field at z = 0, E{z) is the field 

at a distance z into the atomic medium, Im[xxot a i] * s ^ ne 
imaginary part of the linear susceptibility and OD(z) = 
2irN\jli\ 2 1 ki I z/h is the optical depth at z, where ki is 
the wave vector of the field Ei and fli the dipole matrix 
element of the transition. 

We find an OD=20 for light at 780 nm by comparing 
the observed absorption profile in the experiment with 
the theoretical Doppler-averaged absorption profile from 
Eq. © taking into account contributions from the other 
hyperfine levels. There is asymmetry in the absorption 
profile around the detuning for the hyperfine level con- 
sidered in the FWM process 5P 3 / 2 ,P = 4. From the 
absorption data and the model we infer an experimental 
OD for light at 795 nm of 9.8. 

Figure [7] shows the generated light from FWM in 
steady state as a function of the detunings A780 and 
A1324 when the 795 laser is on resonance. Note the asym- 
metry in the intensity of the FWM light from the atomic 
resonances. 

This asymmetry in the resonant structure shows that 
the Rabi frequencies entering into the FWM process are 
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FIG. 7: (Color online) Generated light from FWM in steady 
state as a function of detunings. a) Experimental FWM light 
as a function of 780 and 1324 detunings when the laser at 
795 is locked on resonance, b) FWM light from the model 
taking into account Doppler broadening and propagation ef- 
fects of the fields for incident Rabi frequencies (fli, Q2, Q3) = 
(45, 7, 7) MHz. Color bar in arbitrary units. Note the asym- 
metry in the resonant structure in both cases. 



not symmetric with respect to zero detuning of the lasers. 
This effect modifies the transient behavior of the light 
from FWM at different detunings. The asymmetric be- 
havior of the Rabi frequencies comes from the absorptive 
and dispersive properties of the medium due to hypcrfinc 
levels in the atomic structure of Rb. The model (Fig. [7J 
b) taking into account Doppler effects and absorption 
shows a similar asymmetry in the resonant structure of 
the light from FWM. 



V. RESULTS 

The transient light from FWM in the experiment for 
different detunings shows a shortening of the decay time 



as theoretically predicted in Fig. [6j This is due to 
induced atomic polarization interference in the atomic 
medium [36j . in our particular case from different ve- 
locity classes of atoms in the Doppler-broadened atomic 
ensemble as predicted by Zuo et al. [37]. Different ve- 
locity classes experience different detunings and initial 
phases given by their steady state conditions. The to- 
tal generated FWM field is the sum of all the generated 
fields from different velocity classes. The FWM intensity, 
proportional to the square of the total field and thus to 
the square of the induced atomic polarization, contains 
interference terms from different velocity classes. 

Figure [5] shows a calculation of the intensity of the 
light from FWM with detunings (A795, A1324, A7 80 ) = 
(-60,270,30) MHz and Rabi frequencies (Qi,0 2 ,fi 3 ) = 
(45, 7, 7) MHz for a) the zero-velocity class, b) the sum 
of Doppler-broadened intensities from different velocity 
classes neglecting the induced atomic polarization inter- 
ference in the atomic medium using J \p mn (v)\ 2 P(v)dv 
and c) Doppler broadening including the induced atomic 
polarization interference calculated from the square of 
the Doppler-broadened atomic response using Eq. (|6]) 
1/ Pmn(v)P(v)dv | 2 The Rabi frequencies are consis- 
tent with those estimated from the FWM in the steady 
state regime (see Fig. [7]). 

The transient signal ignoring induced atomic polariza- 
tion interference in the atomic medium in Fig. [8]d has a 
decay time, extracted from a fit of the calculation to a 
simple exponential, of 22.31 ± 3.8 ns. The transient sig- 
nal containing induced atomic polarization interference 
in Fig. [Sfc shows a decay time of 5.82 ± 1.09 ns with a 
modulation frequency of 160 ± 5.4 MHz, both numbers 
extracted from a fit to a simple exponential and from 
the peak frequency in a Fast Fourier Transform (FFT). 
This is consistent with the experimental decay time of 
i~Exp — 6.3 ± 0.09 ns for the experimental Rabi frequen- 
cies shown in Fig. \5jp. The shortening of the decay time 
of the transients with Doppler broadening corresponds to 
a broadening of the FWM spectrum. 

The oscillations in the generated light are the result 
of the average of the oscillatory transients from different 
velocity classes of atoms. The generation of light above 
the steady state level (FWM spike) is due to the total 
contribution to the FWM light from different velocity 
classes whose oscillations phase after the 1324 nm light is 
turned off. At later times, after the in-phase interference, 
the oscillations get out of phase as the different induced 
atomic polarizations from different velocity classes evolve 
in time producing fast damping of the transients. The 
model accounting for Doppler broadening qualitatively 
agrees with the experimental transients shown in Fig. [5] 
for specific detunings of the external excitations. 

Figure [9] shows examples of the modeled FWM tran- 
sients taking into account Doppler broadening using Eq. 
(O and pump depletion using Eq. ©, together with the 
FFT showing their frequency components. 

The oscillation frequency depends on the detunings 
of the lasers. For detunings (A7g 5 , A1324, A7 S0 ) = 
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FIG. 8: (Color online) Theoretical calculation for Transient 
FWM for (A 795 , A1324, A 78 o) = (-60, 270, 30) MHz and Rabi 
frequencies (fh, 02,^3) = (45,7,7). a) Zero-velocity class, 
b) Intensity contributions from all velocity classes excluding 
interference terms, c) Doppler average including interference 
terms. 



(0,270,30) MHz the transient FWM shows an oscil- 
lation frequency at 250 MHz with a half-width at 
half-maximum (HWHM) of 25 MHz. For detunings 
(A 795 , A1324, A 780 ) = (80,2 70,30) MHz the transient 
FWM shows an oscillation frequency at 295 MHz with 
HWHM=30 MHz. The generated light from FWM 
contains contributions from different velocity classes of 
atoms resulting in an averaged oscillation frequency. This 
dephasing mechanism produces broadening as seen in 
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FIG. 9: (Color online) Theoretical calculations of FWM tran- 
sients Doppler-averaged from Eqs. (J3J) and © including 
pump depletion from Eq. ([9]), as a function of detuning A 795 
for (A1324, A 78 o) = (270,30) MHz, together with their Fast 
Fourier Transforms (FFT). a) FWM transients for 795 nm 
resonant and non-resonant light in solid (blue) and dashed 
(red) lines, respectively, b) FFTs with frequency components 
at 250 MHz (HWHM=25 MHz) for resonant light and 295 
MHz (HWHM=30 MHz) for non-resonant light in solid (blue) 
and dashed (red) lines, respectively. The Rabi frequencies 
(fii; £^2 ; ^3) = (45; 7; 7) MHz are consistent with the experi- 
mental parameters. 



Fig. Mb and predicted in Ref. [37) . 

Figure [TU] shows measured examples of the light from 
FWM in the transient regime and the corresponding 
FFT. When the laser detunings are (A795, A1324, A780) = 
(0, 270, 30) MHz the transient FWM contains oscillations 
with a frequency of 185 MHz (HWHM=15 MHz), and for 
the detunings (A 795 , A 1324 , A 780 ) = (80,2 70,30) MHz 
this oscillation changes to 100 MHz (HWHM=20 MHz) 
as shown in Fig. [TUlb. 

Figure [TT] shows the oscillation frequencies /3j from the 
decoupled model of Eq. ([3]) as a function of different de- 
tunings and Rabi frequencies. The continuous lines are 
the results of the model and the points with error bars 
are measurements from FFTs under appropriate experi- 
mental conditions. The horizontal error bars come from 
uncertainty in the frequencies of the lasers. The lasers 
at 780 and 795 nm are locked to atomic absorption lines 
in Rb with linewidths of 6 MHz and frequency-shifted 
using optical modulators (see Sec. We estimate the 
uncertainty of the frequency of the 1324 nm laser from 
the 30 MHz resolution of the wavemeter. 
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FIG. 10: (Color online) Experimental FWM transients as a 
function of detuning A 795 for (A1324, A 78 o) = (270,30) MHz 
and FFTs for a situation similar to Fig. [21 a) FWM transients 
for 795 nm resonant and non-resonant light with data points 
joined by solid (blue) and dashed (red) lines, respectively, b) 
FFTs with frequency components at 185 MHz (HWHM=15 
MHz) for resonant light and 100 MHz (HWHM=20 MHz) for 
non-resonant light with data points joined by solid (blue) and 
dashed (red) lines, respectively. 



This simplified model shows some characteristics of the 
generated light. The transients should oscillate at the 
natural frequencies of the induced atomic polarization in 
the medium, however, we observe some discrepancy be- 
tween the experimental observations and the expected 
frequencies from the model. These discrepancies could 
come from the fact that this simplified model does not 
contain information about pump absorption or transverse 
effects. Pump depletion could produce abrupt changes in 
the Rabi frequencies entering into the FWM process as 
the lasers are scanned across atomic resonances because 
of the presence of other hyperfine levels. Transverse ef- 
fects in the atomic medium could modify the wave vectors 
of the pumps as a function of frequency and change the 
phase matching condition required in FWM 

The asymmetry with respect to the 795 and 780 de- 
tunings in the experimental oscillation frequencies might 
be related to the asymmetric behavior of the observed 
FWM light in the steady state in Fig. [7k, The asym- 
metric absorption of the pump beams around the partic- 
ular hyperfine levels predicts an asymmetry in the FWM 
in the steady state regime (see Fig. [7p). However, a 
model fully describing the asymmetric behavior of the 
transient oscillations in Fig. [TT] might require solving the 
coupled density matrix and field propagation equations 
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FIG. 11: (Color online) Experimental (points with error 
bars) and modeled (continuous lines) oscillation frequencies 
Pjk (j 7^ fc) as a function of detunings. a) As a function 
of A 780 with (A 795 ,Ai324) = (0,270) MHz. b) As a func- 
tion of A 795 with (A 78 o,Ai324) = (30,270) MHz. c) As a 
function of A1324 with (A 780 ,A 795 ) = (30,0) MHz. The 
model ignores motion of the atoms. The Rabi frequencies 
(^1,^2,^3) = (95,20,20) MHz. The different colors of the 
lines correspond to oscillation frequencies from specific pairs 
of quasienergies in Eq. ([5]) . 



to account for possible coupling between the pump fields 
in the medium producing focusing and defocusing as well 
as nonlinear refractive index effects modifying the phase 
matching for the fields 

The location of the oscillations (data points) shows 
that the system changes among the different oscillation 
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frequencies available. The data comes from the FFT of 
the time response and one oscillation dominates in every 
case. Further investigation will be necessary to under- 
stand this effect in the FWM process. 

VI. CONCLUSIONS 

Coherences, such as those created in the FWM pro- 
cess, are delicate. This paper shows the atomic polar- 
ization coherences that contribute to the FWM and are 
visible through the time response of the system. They 
are responsible for transients, but their fragility is evi- 
dent in the qualitative and quantitative behaviors that 
strongly depend on detunings and intensities. Our work 
also shows that the presence of absorption and disper- 
sion with the Doppler-broadened medium significantly 
changes the coherences and their dynamics: oscillations 
and decay time. 

The observations of pure decay agree with the model 
for resonant excitation and low power; changes in the 
detuning bring oscillations and jumps above the steady 
state value of the FWM signal which come from induced 
atomic polarization interference in the atomic medium. 



The total transients result from contributions of differ- 
ent velocity classes of atoms in the Doppler-broadened 
medium; this produces induced atomic polarization in- 
terference generated from different velocities broadening 
the FWM spectrum. The model predictions for the am- 
plitudes of the transients in FWM are in qualitative agre- 
ment with the experimental observations. The frequency 
structure, related to the quasienergies of the system, is 
closer quantitatively to the model. 

This work provides a better understanding of coher- 
ent effects and energy exchange processes in multilevel 
atomic systems and how they are affected due to broad- 
ening as a source of decoherence. We find that, in this 
Doppler-broadened system, it is possible to obtain co- 
herent processes with decay times of the order of the 
radiative lifetimes for transitions in the diamond config- 
uration. 
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